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Abstract. Recently it was shown that the exact cosmological solutions known 
as the vacuum plane-wave solutions are late-time attractors for an open set 
of the spatially homogeneous Bianchi universes containing a non-inflationary 
7-law perfect fluid. In this paper we study inhomogeneous perturbations of 
these plane-wave spacetimes. By using expansion-normalised scale-invariant 
variables we show that these solutions are unstable to generic inhomogeneous 
perturbations. The crucial observation for establishing this result is a diver- 
gence of the expansion-normalised frame variables which ultimately leads to 
OO . unstable modes. 
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O . 1- Introduction 

iy-) . Since the advent of general relativity researchers have sought physical relevant 

solutions to the Einstein field equations. However, due to the complicated nature 
of the field equations symmetries are usually imposed in order to make the field 
equations more tractable; for example, the spatially homogeneous Bianchi models 



o 



in (3 + l)-dimcnsional cosmology are of particular interest |2 |3J 0]. One of 



the most fascinating and interesting solutions to these field equations are the so- 



5-H . called vacuum plane-waves. These solutions have some extraordinary physical and 



mathematical properties which have drawn the attention of a broad spectrum of 
researchers. For the Bianchi models these plane-waves play a particular role; they 
are late-time attractors for non-inflationary non-compact perfect fluid cosmologies 
of Bianchi class B. 



& ' The spatially homogeneous (SH) Bianchi class B models contain types IV, V, 

VI/,, and VII/j. Bianchi models are referred to as orthogonal cosmologies when 
the fluid flow is orthogonal to the group orbits. Otherwise they represent 'tilted' 
models. The types VI/j and VII^ are of non-zero measure in the space of all spa- 
tially homogeneous models and are thus what might be called 'general' type. The 
type Vllh models are of particular physical interest because they contain open 
Friedmann- Robertson- Walker (FRW) models. It is known that in general relativity 
all orthogonal class B perfect fluid models, with equation of state parameter 7 < 2, 
expand indefinitely into the future and an open set of models are asymptotic in 
the future to a vacuum plane wave (self-similar) solution E>| . In the stiff fluid 
case, apart from a set of measure zero, the asymptotic form of the general type 
VII/! stiff perfect fluid solution corresponds to the vacuum type VII^ plane wave 
solution first derived by Doroshkevich, Lukash and Novikov [Jf . Orthogonal mod- 
els are special; it is therefore important to consider the effect that tilt has on the 
asymptotic behaviour of the class B models. 
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Until recently, only partial results regarding tilted Bianchi models were known 
[HUH]. However, a detailed analysis of the tilted Bianchi models shows that many 
of the vacuum plane waves are stable with regards to more general tilted matter 
|101 ITT] . This analysis includes keeping track of the asymptotic tilt velocity of the 
perfect fluid which, in general, can have different kinds of complicated asymptotic 
behaviour; e.g., attracting points, closed orbits and even attracting tori. Notwith- 
standing this complex behaviour of the tilt velocity, the plane-wave geometry are 
asymptotically stable for types IV, VII^, and an open set of type VI^ models. Since 
the types VII/j and VI/j are among the most general Bianchi models, we can draw a 
remarkable conclusion: In the class of all spatially non-compact Bianchi universes 
with a non-inflationary "/-law perfect fluid (not necessary co-moving) , there exists 
an open set of universes asymptoting towards the vacuum plane-waves at late times. 

As stressed in 0] , an important mathematical link exists between the various 
classes in the state space hierarchy. The physical state of a cosmological model at 
an instant of time is represented by a point in the state space, which is finite di- 
mensional for SH models and infinite dimensional otherwise. This structure opens 
the possibility that the evolution of a model in one class may be approximated, 
over some time interval, by a model in a more special class. Thus it is plausible 
that understanding the dynamics at one level of complexity, such as for example 
the spatially homogeneous Bianchi type VI^ and VII^ models, will shed light on 
the possible dynamical behaviour at a higher level, such as special classes of inho- 
mogeneous cosmological models. 

To investigate this further, in this paper we will consider inhomogeneous per- 
turbations of the Bianchi type VII;, plane-wave spacetimes. In particular, we will 
determine whether the plane- wave solutions are future attractors for models of non- 
zero measure in general relativity. However, our analysis indicates that this is not 
so; the models are unstable with respect to inhomogeneous perturbations. The crit- 
ical unstable modes we find correspond to inhomogeneous perturbations orthogonal 
to the propagation of the gravitational wave. For these plane-wave spacetimes the 
direction along the propagation of the wave expands linearly in cosmological time, 
while the orthogonal directions expands more slower. The instability of these mod- 
els stems from this anisotropic expansion; the slower expanding directions do not 
expand rapidly enough to prevent the forming of inhomogeneities. 

In order to understand this instability, it is illustrative to compare the expansion 
with a flat FRW model. The flat FRW model expands as a(t) oc t p , where t is the 
cosmological time. Inflation occurs for p > 1, while p < 1 leads to instability with 
regards to inhomogeneous perturbations jo]. For the plane- wave spacetimes, one 
direction expands with p = 1 (the threshold value), while the others expand with 
p < 1 . Hence, we expect instability in two of the directions (the third is a marginal 
case). 



2. Inhomogeneous cosmologies 

Recently, a formalism for studying the dynamics of inhomogeneous cosmologies 
(Go models) by employing expansion-normalised scale-invariant variables was de- 
veloped |12| . This formalism uses orthonormal frames to write the Einstein field 
equations as an autonomous system of evolution equations and constraints. 
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The orthonormal frame {eo,e a } can be expressed in local coordinates as 

(1) e = jj(dt-N i d i ), e a = eX 
We introduce the Hubble-normalised frame as follows: 

(2) d = -^e , d a = —e a . 

It is also necessary to introduce the deceleration parameter q and the spatial Hubble 
gradient r Q : 

(3) q + l = -^d H, r a = -^d a H. 

For the purposes of this paper, it is useful to introduce the separable volume gauge: 

(4) N = H-\ AP = 0, U a = r a . 
This enables us to write 

(5) d =d t , d a =E a %, 
where = e*/H. 

The Hubble-normalised state vector for Go cosmologies with a 7-law perfect fluid 
can then be given by 

(6) X= [E a \ r a , £ Q/3 , A a , N af3 , ft, v a f , 

where Yi a p is the expansion-normalised shear, A a and N a p are the connection vari- 
ables, and f2 and v a are the energy density and the tilt of the matter, respectively. 
The evolution equations and the constraints can now be written down in terms of 
X (see 521)- The equations can be written in the form 

(7) a t x = F(x,9 l x,a,a J x), 

(8) = C(X,diX), 

where eqs.0 are the evolution equations and eqs. © are the constraints. An 
interesting observation is that apart from the evolution equations for r a , the spatial 
derivatives appear linearly. Moreover, the only equations involving second-order 
spatial derivatives are those of the Hubble gradient, r a , which contain terms of the 
form didjr a . 

2.1. The plane waves. In our case we will be interested in a particular set of 
exact solutions, namely the vacuum plane-wave solutions. In general the plane- 
wave solutions have a 2-dimensional Abelian group acting freely on the spatial 
hypersurfaces 1 ; i.e., they are special exact G2 solutions to the equation of motion. 
In particular, by choosing the vector r a appropriately, the G2 models have r a — 
(r i,0,0) (see for details). The exact G2 plane- waves solutions can easily be 
found; however, we will only consider a special class of them here. The class we 
are considering is the spatially homogeneous plane-wave solutions of Bianchi type 
Vllfc. These plane- wave solutions are given by 

ds 2 = e 2T (-dr 2 + dx 2 ) + e 2s(x+T) [e^cosudy - sinvdz) 2 + e _/3 (sint)d?/ + cosvdz) 2 ] 

where v = b(x + t) and 6, j3 and s are constants satisfying b 2 sinh 2 f3 = s(l — s), 
< s < 1. The Killing vectors d y and d z span the wave-fronts and are the Killing 

^Cosmological models which have a 2-dimensional Abelian group acting freely on the spatial 
hypersurfaces are usually referred to as G2-models pj- 
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vectors that survive the generalisation to the G 2 models. We will see later that the 
critical unstable modes correspond to inhomogeneous perturbations in the (y, z)- 
plane so we do not expect the specialisation to spatially plane-wave homogeneous 
models to be critical for our result. 

3. Perturbing the plane wave 

The spatially homogeneous plane waves of Bianchi type VII^ are given in the 
expansion- normalised variables by (the notation used is the same as in |l()l 111) ) 

(9) q = 2 X , £+ = -X, A a = ASi, A=l- Xl Z 2 _ = 7V| 3 = x (l - x ), 
N = XN 23 , where < x < 1, A 2 > 1, (1 - X ) = 3/i(A 2 - l) X 

£ 23 - Eia = Eia = iV_ = Nu = N 12 = N 13 = r a = Q = 0. 

In the following we will assume that O = (vacuum) 2 . Moreover, we will also chose 
the following gauge: 

(10) Ri = V3A£_, R 2 = -V3Si3, #3 = V3E 12 . 

An alternative choice of gauge could be, for example, one where we eliminate some 
of the expansion-normalised variables; however, the above choice is useful because 
it simplifies the equations of motion significantly. We will also assume that the 
Codazzi equations has been used to eliminate the variables £12 and £13. 

3.1. Showing instability. Let us first establish the fact that these plane waves 
are unstable to inhomogeneous perturbations. Consider the variable Nu, which 
vanishes for the background plane-wave spacetime. For the above choice of gauge 
the evolution equation for N±i is 

(11) atJVii = («-4E+)JVn. 

At the linearised level, we get dtNn — 6 X Nn, and hence this variable is unstable. 
So for models that allow for a non-zero Nu, the plane waves are unstable. The 
constraints governing the N a p variables, are the Jacobi constraint equations: 

(12) = (0/j - 173) (N aP + e a ^A^) - 2A p N aP . 

It is clear that the general inhomogeneous models allow for a non-zero N\\\ hence, 
the plane-wave spacetime is unstable to generic inhomogeneous perturbations. 
The component Nu is the 1-componcnt of the commutator between d 2 and Q 3 : 

(13) tfn = ([02,0a])!, 

which means that a pertubation in Nu stems from a spatial perturbation of the 
vectors d 2 and 03. Thus to understand where this instability comes from, it is of 
importance to study the differential operators d 2 and 3 , and the Hubble gradients 
r a = —(l/H)d a H. Note that in the Abelian case (i.e., the G 2 models) Nu = 0. 



The assumption Q = is not crucial for our analysis since we find that the models are 
unstable. Including matter will only bring more degrees of freedom, and hence more perturbative 
modes, into the analysis. 
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3.2. Interpretation. The frame variables, EJ, decouple for the spatially homo- 
geneous models where eo is orthogonal to the type VII^ hypersurfaces. For the 
background type VII^ plane waves, the frame variables can be taken to have the 
form 

Ej = (1,0,0) 

E 2 < - (o,E 2 y,E 2 >) 

(14) Ej = {0,E*,E 3 *), 

where the matrix (£J a l ) is invertible. Introducing the complex frame variable 

Z* = Ej + 1 (l + iy/V - l) Ej, i = y,z, 
the evolution equations and the commutator relations, 

(15) d t Ej = (qS^-X^ + e^R^Ej 

(16) = 2(d la -r [a -A la )Efi-e afiS N*iE 7 
imply that for the plane waves 

d t Z l = {2> X -iuj)Z\ 

(17) d x Z l = (A-iuj)Z\ 



where ui = y / 3x(l — x)(^ 2 — !)• Hence, cx exp[3x£+Ar— iu(t+a;)]. This implies 
that the frame variables diverge for the plane waves in question. Note that for a SH 
cosmology, we have that <9iX = 0, and thus d a ~K = 0. So in these models the frame 
variables decouple from the equations of motion. However, as we perturb these 
models inhomogeneously, e? Q X is no longer necessarily zero and the divergence of 
the frame variables may become important. The above shows that the differential 
operators &2 and c?3 carry a divergent term e 3x * which, in turn, indicates that the 
plane waves are unstable with regard to inhomogeneous perturbations in the d 2 
and c?3 directions (orthogonal to the propagation of the wave). 

In the following analysis we will choose Z y — exp[3%t + Ax — iu>(t + x)] = Z, 
and Z z = iZ. It is also advantageous to introduce the following complex variable, 

r = r 2 + j(l + i^J\ 2 - l) r 3 

along with the operators 

d = Z*d l= Z(dy+ld Z ), 

(18) B = z l d t = z(d v -id z ). 

To investigate the (in)stability of these spacetimes with respect to general in- 
homogeneous perturbations we will linearise the equations of motion around the 
exact plane-wave solutions (assuming that terms like (X — X ), c? Q X, d a dpX., etc. 
are small). The interesting modes stem from the Hubble gradients r\ and r, so 
to understand what is happening we will linearise these equations explicitly. The 
evolution equation for r a and the vorticity constraint read 

(19) 8 t r a = (q5^ a -^ a +E^If)rp + d a q 

(20) = [e afh (dp - A p ) - N a ~<] r 7 . 
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The deceleration parameter, q, for vacuum is given by 

q = 2Y? + l i A a r a -^d a r a . 

From the linearised vorticity constraint, dr = dr, and using [d, d] — 0, we can 
write 

(21) dq = 2dY? + \{2Adn - dd x n) - ^(d 2 2 + tifrr. 
The vorticity constraint also gives the following linear relation: 

(22) dn = (d x - A)r + M r r + M 7 r, 

where M r and M 7 are constants. However, for simplicity, we will leave this constraint 
unsolved and rather impose it on the initial conditions. 

At this stage it is useful to introduce an orthonormal set of functions, / k (y, z), on 
the flat torus spanned by the vectors d y and d z . These functions can be chosen to 
be proportional to the usual Fourier modes e ~ l ( kyV+kzZ \ and hence, S/ k — —ikZfk, 
where k = k y + ik z . Defining 

(23) « 2 = E 2 % = -^=(kZ-kZ), 



kZ 



(24) k 3 = E 3 % = 7^==j [(l - ^A 2 - l) kZ - (l + iV\ 2 - 1 

and k 2 = k\ + k§ > 0, we obtain 

(25) (dl + dj)h = -« 2 / k . 

Here, n 2 is a quadratic monomial in Z and Z, and therefore diverges as k 2 oc e 6x *. 
By expanding all the linearised variables as 

(26) X = ^X k (i,aO/ k (y,z), 



4(£ Q/3 ) Ef + H(i _ X ) ri k _ l dx r 1M 



we obtain 

<9 t r k = ^3% - «w + ^K 2 ^j r k - zkZ 
(27) 

We note that we can write 

(28) (|r k | 2 )' - (& X + 1^ |r k | 2 - 2Rc(*kZp k r k ), 

where p k is a linear combination of Hff n ik and d x r\^. For k ^ 0, there exist 
{^minW.^L^)} > such that 

e 6xt K 2 nin (a;) < k 2 < e 6x *K 2 nax (x). 

Let us also introduce 

(29) I = | / K 2 dt, 



to 



which obeys the bounds 



K min( X ) I 6 X t _ e 6xta \ < J < K max( x ) / 6 X t _ e §Xto\ 
9X " 9 X ^ ' 
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The equation for n has a similar k 2 term, so taking into account the vorticity 
constraint, we would naively expect that |rk| 2 oc e At+/ for a constant A, and hence, 
rk is diverging. 

We are thus led to conclude that the plane waves are unstable with respect to 
inhomogeneous perturbations of the Hubble gradient r. 

4. Discussion 

In this paper we have shown that the plane-wave spacetimcs of Bianchi type 
Vllh are unstable with respect to inhomogeneous perturbations. The cause of this 
instability ultimately stems from the behaviour of the frame variables E^, which 
are unbounded into the future for the plane-wave solutions. For the spatially ho- 
mogeneous solutions these frame variables decouple from the equations of motion; 
however, in the general Go models this unboundedness of the frame variables is cru- 
cial. As we pointed out earlier, the expansion-normalised differential operators Q a 
carry this divergence and are thus in some sense intrinsically unstable. In partic- 
ular, we showed that the Hubble gradient, r a = —(l/H)d a H, has two unbounded 
components corresponding to directions orthogonal to the propagation of the wave. 

We also point out that the limit \ — > corresponds to the vacuum Milne model. 
The above analysis has excluded this case since the eigenvalue (3% — iu) gets a 
troublesome zero real part. In this case a separate analysis is required. Pertubations 
of the Milne universe have been studied in which shows that this model is also 
unstable. 

For the G2 models, r a = (n, 0, 0) and hence the critical unstable modes obtained 
here are excluded. A non-zero n is still allowed, and from the equation of motion 
we note that there is a 'zero eigenvalue' corresponding to the evolution equation 
for r% . This is really no surprise since we know that there exist more general plane 
wave solutions with n 7^ 0. It is not known what role these G2 plane waves play 
in the evolution of general G2 cosmologies. 

Another point worth mentioning is the choice of temporal gauge. In this work 
we used the separable volume gauge, defined in cq.(0J; however, another useful 
gauge is the synchonous gauge. Also in this gauge we can show that the plane 
waves are unstable; in particular, eo.^lljl remains on the same form. This shows 
the robustness of our result with regards to the temporal gauge. 

Recently, Barrow and Tsagas have studied linearised perturbations of the 
same plane-wave spacetimes. However, they considered the subcase in which the 
Weyl tensor does not vary along the directions d y and d z . This assumption alone 
effectively reduces the model to a Gi model and thus, as we pointed out eariler, 
the unstable critical mode we found here is not included. 

In the analysis the unboundedness of the frame variables E^ was critical for the 
unstable nature of these models. We may thus ask ourselves: Will all models which 
have unbounded frame variables E£ be unstable against generic inhomogeneous 
perturbations? Future research may shed some light on this issue. 
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